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A truncation scheme for the Dyson-Schwinger equations of QCD in Landau gauge 
is presented which implements the Slavnov-Taylor identities for the 3-point vertex 
functions. Neglecting contributions from 4-point correlations such as the 4-gluon 
vertex function and irreducible scattering kernels, a closed system of equations for 
the propagators is obtained. For the pure gauge theory without quarks this system 
of equations for the propagators of gluons and ghosts is solved in an approximation 
which allows for an analytic discussion of its solutions in the infrared: The gluon 
propagator is shown to vanish for small spacelike momenta whereas the ghost 
propagator is found to be infrared enhanced. The running coupling of the non- 
perturbative subtraction scheme approaches an infrared stable fixed point at a 
critical value of the coupling, a c ~ 9.5. The results for the propagators obtained 
here compare favorably with recent lattice calculations. 



1 Introduction 

The main objective of this talk is to present results for the basic QCD propa- 
gators and for the running coupling. Knowledge of these quantities may yield 
further insight into the physical nature of confinement. One appealing possi- 
bility to describe confinement is to ascribe the suppression of the emission of 
colored states from color-singlet states to infrared divergencies in the elemen- 
tary correlators of QCD. 

Thus, to study the infrared behavior of QCD amplitudes non-perturbative 
methods are required, and, since divergences are anticipated, a formulation in 
the continuum is desirable. Both of these are provided by studies of trun- 
cated systems of Dyson-Schwinger equations (DSEs) , the equations of motion 
of QCD Green's functions. Typically, for their truncation, additional sources 
of information like the Slavnov-Taylor identities, entailed by gauge invariance, 
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are used to express vertex functions in terms of the elementary 2-point func- 
tions, i.e., the quark, ghost and gluon propagators. Those propagators can 
then be obtained as selfconsistent solutions to non-linear integral equations 
representing a closed set of truncated DSEs. Some systematic control over the 
truncating assumptions can be obtained by successively including higher n— 
point functions in selfconsistent calculations, and by assessing their influence 
on lower n-point functions in this way. Until recently all solutions to truncated 
DSEs of QCD in Landau gauge n eveii in absence of quarks, relyed on neglecting 
ghost contributions completelyEl'B'Hu. While this particular problem is avoided 
in ghost free gauges such as the axial gauge, in studies of the gluon DSE in this 
gauge Q, the possible occurrence of an independent second term in the tensor 
structure of the gluon propagator has so far been disregardedla. In fact, if the 
complete tensor structure of the gluon propagator in axial gauge is taken into 
account, one arrives at equations of no less complexitiy than the ghost-gluon 
system in the Landau gauge. 

In this talk a simultaneous solution of a truncated set of P.SEs for the prop- 
agators of gluons and ghosts in Landau gauge is presented LrEl. An extension 
to this selfconsistent framework to include quarks is possible and subject to 
further studies. The behavior of the solutions in the infrared, implying the ex- 
istence of a fixed point at a critical coupling a c ~ 9.5, is obtained analytically. 
The gluon propagator is shown to vanish for small spacelike momenta in the 
present truncaticu^scherne. This behavior, though in contradiction with previ- 
ous DSE studiestloEKa, can be partially understood from the observation that, 
in our present calculation, the previously neglected ghost propagator assumes 
an infrared enhancement similar to what was then obtained for the gluon. 

In the meantime the qualitative behavior obtained in the studies to be 
reported here is supported by investigations of the coupled gluon ghost DSEs 
using bare vertices with □ and without E3 angle approximation. As expected 
however, the details of the result depend on the approximations employed. 

2 The set of truncated gluon and ghost DSEs 

Besides all elementary 2-point functions, i.e., the quark, ghost and gluon prop- 
agators, the DSE for the gluon propagator also involves the 3- and 4-point 
vertex functions which obey their own DSEs. These equations involve suc- 
cessively higher n-point functions. A first step towards a truncation of the 
gluon equation is to neglect all terms with 4-gluon vertices. These are the 
momentum independent tadpole term, an irrelevant constant which vanishes 
perturbatively in Landau gauge, and explicit 2-loop contributions to the gluon 
DSE. The latter are subdominant in the ultraviolet and will thus not affect the 
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behavior of the solutions for asymptotically high momenta. In the infrared it 
has been argued that the singularity structure of the 2-loop terms does not 
interfere with the one-loop terms Ell. Without contributions from 4-gluon ver- 
tices (and quarks) the renormalized equation for the inverse gluon propagator 
in Euclidean momentum space is given by 

D~i{k) = Z 3 D^{k) +g 2 N c Z^ J J^jL 

x r M PQ ( fc ; -P> ~i) D a p{q)Dpa{p) I> CT „(g,p, -k) 

/i4- 
—^i qil D G (p)DG{q)G v {-q,p), (1) 

where we use positive definite metric, g^ v — 5^. Color indices are suppressed 
and the number of colors is fixed, N c = 3. Furthermore p — k — q, D a and P 1 
are the tree level propagator and 3-gluon vertex, Dq is the ghost propagator 
and T and G are the fully dressed 3-point vertex functions. The equation for 
the ghost propagator in Landau gauge QCD, without any truncations, is given 
by 

D G \k) = -Z 3 k 2 + g 2 N c Z 1 (2) 

77^ v4 ik v D^ u (k - q) G v (k, q) D G {q) . 

The renormalized propagators for ghosts and gluons and the renormalized cou- 
pling are defined from the respective bare quantities by introducing multiplica- 
tive renormalization constants, Z^Dq ■— D G , Z^D^ := and Z g g := go. 

3 /2 ~ 1/2 — — 

Furthermore, Z\ = Z g Z^ , Z\ = Z g Z^ Z3, and we use that Z\~\ in Landau 
gauge. The ghost and gluon propagators are parameterized by their respective 
renormalization functions G and Z, 



n •/• G ^ n m ( A ^k v \Z(k 2 ) 
D c(k) = p-, £>„*(*) = I V--p-)-p-- ( 3 ) 



In order to arrive at a closed set of equations for the functions G and Z, 
we use a form for the ghost-gluon vertex which is based on a construction 
from its Slavnov-Taylor identity (STI) which can be derived from the usual 
Becchi-Rouet-Stora invariance neglecting irreducible 4-ghost correlations in 
agreement with the present level of truncation a. This together with the sym- 
metry of the ghost-gluon vertex fully determines its form at the present level 
of truncation. There are no undetermined transverse terms in this case: 

~2\ 



M = ^^ + ^(^p)-i| (4) 



3 



With this result, we can constmct the 3-gluon vertex according to procedures 
developed and used previously 113, 

IW(p, q, k) = ^A + (p 2 , q 2 ; k 2 ) 5^ i(p - q) p (5) 

+ U_ (p 2 ,q 2 ; k 2 ) 8» v i(p + q) p + A ~^ff) 
L p q 

x [Sftvpq - p v q p ) i(p - q) p + cyclic permutations , 

JUtf a 2. h 2, G(k 2 )G(g 2 ) G{k 2 )G{p 2 ) 
A±[P > q '* >- G(p 2 )Z(p 2 ) ± G(q 2 )Z(q 2 ) ■ 

Some additionally possible terms, transverse with respect to all three gluon 
momenta, cannot be constrained by its STI and are thus disregarded. 

We solve the coupled system of integral equations of the present trunca- 
tion scheme in an one-dimensional approximation: For integration momenta 
q 2 > k 2 we use the perturbative anticipation that the functions Z and G are 
slowly varying with their arguments, and that we may thus replace all argu- 
ments by the integration momentum q 2 . This assumption ensures the correct 
leading perturbative behavior at one-loop_ level An analogous assumption 
underlies the Mandelstam approximation EJH'B'Q. Infrared enhanced solutions 
tend to invalidate this assumption, in particular for small q 2 < k 2 , however. 
For integration momenta q 2 < k 2 we therefore use the angle approximation 
replacing G{{k — q) 2 ) — > G{k 2 ) and similarly for Z , which preserves the limit 
of the integrands for q 2 — > 0. The DSEs (Q) and (||) then simplify to 

1 7 ^ 7 9 2 j f^dq 2 (7q* 17 q 2 9\ 

— — = Z 3 + Zt—ri / — -— - — — - - )Z(q )G{q ) 



d 4%G{ki)G{?) - W) + i (6) 



16tt 2 4 



(7) 



We introduced an 0(4)-invariant momentum cutoff A uv to account for log- 
arithmic ultraviolet^ divergences which are absorbed by the renormalization 
constants Z3 and Z3. The details of the renormalization and the numerical 
procedure are given elsewhere El. 
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The angle approximation as we use it for small integration momenta was 
also used in a very .recent study of the coupled system of ghost and gluon DSEs 
with bare vertices B. However, using this approximation for arbitrary q 2 (i.e. 
also for q 2 > k 2 ) one does not recover the renormalization group improved 
one-loop results for asymptotically large momentaEl. 

For a reduced set of equations the effect of the angle approximation has 
been demonstrated to be only quantitative, the qualitative infrared behavior of 
the solution remaining the samellj. It will nevertheless be important to assess 
the sensitivity of the results to the modified angle approximation further in 
future. 

To deduce the infrared behavior of the propagators here, we make the 
Ansatz that for x := k 2 — > the product Z(x)G(x) — * cx K with n ^ and 
some constant c. The special case k = leads to a logarithmic singularity in 
Eq. (0) for x —>■ which precludes the possibility of a selfconsistent solution. In 
order to obtain a positive definite function G(x) for positive x from an equally 
positive Z(x), as x — > 0, we obtain the further restriction < K < 2. Eq. (fy 
then yields, 

(/7o G Q-0) V 1 *"" => (8) 
^)-(5 2 7o G Q-0)cV K , (9) 

where 7g = 9/(647r 2 ) is the leading order perturbative coefficient of the anoma- 
lous dimension of the ghost field. Using (||) and (||) in Eq. (^|), we find that 
the 3-gluon loop contributes terms ~ x K to gluon equation for x — >■ while 
the dominant (infrared singular) contribution arises from the ghost -loop, i.e., 

Comparing this to ([)]) we obtain a quadratic equation for n with a unique 
solution for the exponent in < n < 2: 

..etww^M,. (10) 

The leading behavior of the gluon and ghost renormalization functions and 
thus of their propagators is entirely due to ghost contributions. The details 
of the approximations to the 3-gluon loop have no influence on the above 
considerations. Compared to the Mandelstam approximation, in which the 
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Figure 1: The present result compared to the world's experimental data according to ref. 
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3-gluon loop alone determines the infrared .behavior of the gluon propagator 
and the running coupling in Landau gaugeBBoa, this shows the importance of 
ghosts. The result presented here implies an infrared stable fixed point in the 
non-perturbative running coupling of our subtraction scheme, defined by 



a s {s) 



An 



Z(s)G 2 (s) 



16tt 
~9~ 



9.5 



(11) 



for s — > 0. This is qualitatively different from the infrared singular coupling of 
the Mandelstam approximation a. 

The momentum scale in our calculations is fixed from the phenomenologi- 
cal value as(Mz) = 0.118 at the mass of the Z-boson. The ratio of the Z— to 
the r-mass, Mz/M T ~ 51.5, then yields as(M T ) — 0.38. In Fig. |l]we compare 
our result far-the running coupling to the experimental data summarized in 
fig. 6 of ref. £3. The enhancement at lower momenta can be attributed to the 
omission of quarks in our calculation (Nf = 0). 



3 Comparison to lattice results 

It is interesting to compare our solutions to recent lattice results available for 
the gluon propagatortil and for the ghost propagatort3 using lattice versions to 
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Figure 2: The numerical result for the gluon projpagator from Dyson-Schwinger equations 
(solid line) compared to recent lattice calculations uf . 

implement the Landau gauge condition. In Fig. || we compare our solution for 
the gluon propagator to the most recent data for the lattice gluon correlator O 

In Fig. U we compared our infrared enhanced ghost propagator to corre- 
sponding recent results Ej. It is quite amazing to observe that our solution fits 
the lattice data at low momenta significantly hatter than the fit to an infrared 
singular form Do{k 2 ) = c/k 2 + d/k 4 given inE3. We therefore conclude that 
present lattice calculations confirm the existence of an infrared enhanced ghost 
propagator of the form Dq ~ l/(k 2 ) 1+K with < k < 1. This is an interesting 
result for yet another reason: In this calculation the Landau gauge condition 
was supplemented by an algorithm to select gauge field configurations from the 
fundamental modular region which is to avoid Gribov copies. Thus, our results 
suggest that the existence of such copies of gauge configurations might have 
little effect on the solutions to Landau gauge DSEs. This could also explain the 
similarity of our solutions to the infrared behavior obtained by D. ZwanzigerliU 



4 Summary 

In summary, we presented a solution to a truncated set of coupled Dyson- 
Schwinger equations for gluons and ghosts in Landau gauge. The infrared 
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Figure 3: The numerical result for the ghost propagptpr from Dyson— Schwinger equations 
(solid line) compared to lattice data from Fig. 1 in inEj and a fit as given by these authors. 

behavior of this solution, obtained analytically, represents a strongly infrared 
enhanced ghost propagator and an infrared vanishing gluon propagator. 

Using bare vertices qualitatively similar results are obtained in an inves- 
tigation of the coupled gluon and ghost DSEs. Relying on an angle approxk 
mation one obtains k ~ 0.77 for the exponent of the infrared leading-.term □ 
whereas without angle approximation the exponent changes to n = lE2l. 

A highly infrared singular behavior for the ghost propagator in Landau 
gauge has also been suggested from studies of the influence of a complete 
gauge fixingtLj. Our results for the propagators, in. particular f° r the ghost, 
compare favorably with recent lattice calculations oEj. 

The Euclidean gluon correlation function presented here can be shown to 
violate reflection positivityu, which is a necessary and sufficient condition for 
the existence of a Lehmann representation^. We interpret this as representing 
confined gluons. In order to understand how these correlations can give rise to 
confinement of quarks, it will be necessary to include the quark propagator. 

The size of the coupling at the fixed point, a c » 9.5, is however, a good 
indication that dynamical chiral symmetry breaking will be generated in the 
quark DSE. The existence such an infrared fixed point is in qualitative disagree- 
ment with previous-studies of the gluon DSE neglecting ghost contributions in 
Landau gauge aBBa This shows that ghosts are important, in particular, at 
low energy scales relevant to hadronic observables. 
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